function [output] = cheb1(wp, ws, rp, rs, input)
[N,wn]=cheb1ord(wp,ws,rp,rs);
[b,a]=cheby1(N,rp,wn);
output=filter(b,a,input);
